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Abstract 



This paper studies the electroweak production of all possible four- 
fermion states in e + e~ collisions with non-standard triple gauge boson 
couplings. All CP conserving couplings are considered. It is an ex- 
tension of the methods and strategy, which were recently used for the 
Standard Model electroweak production of four-fermion final states. 
Since the fermions are taken to be massless the matrix elements can be 
Q | evaluated efficiently but certain phase space cuts have to be imposed 

to avoid singularities. Experimental cuts are of a similar nature. With 
QJ the help of the constructed event generator a number of illustrative 

^ results is obtained for VF-pair production. These show on one hand 

the distortions of the Standard Model angular distributions caused by 
either off-shell effects or initial state radiation. On the other hand, 
also the modifications of distributions due to anomalous couplings are 
presented, considering either signal diagrams or all diagrams. 
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1 Introduction 



Recently, the electroweak four-fermion production processes relevant for LEP2 
and beyond have been studied in a number of ways. One of the objectives is 
to obtain a description of PU-pair production better than an on-shell treat- 
ment with PU-decay products attached to it. Thus all recent papers contain 
finite width effects. Some papers only include the three diagrams leading to 
PU-pair production, others include all diagrams giving a specific four-fermion 
final state. Most of them include some form of initial state QED radiative 
corrections. There are semi-analytical methods [1, 2] and Monte Carlo ap- 
proaches [3]-[10]. The former can only give distributions in the virtualities 
of the W's, but no fermion distributions. The latter can produce any wanted 
distribution. 

Among the various Monte Carlo treatments we mention in particular 
the program EXCALIBUR, since it aims both at completeness and speed. 
All diagrams for any four-fermion final state are included and a relatively 
fast calculation is achieved by assuming massless fermions and by using a 
multichannel approach to generate the phase space. The details are given 
in [9], whereas the treatment of initial state radiation (ISR) can be found 
in [10]. 

One of the objectives of LEP2 and future electron-positron colliders is a 
test of non-abelian triple gauge boson couplings. A way to quantify devia- 
tions from the Standard Model (SM) Yang-Mills couplings is to set experi- 
mental limits on anomalous couplings. Many discussions of the latter can be 
found in the literature, see e.g. [11, 12, 13]. Theoretical arguments, which 
reduce the a-priori large number of non-standard couplings are discussed 
in [13]. In order to investigate the experimental possibilities to measure lim- 
its on anomalous couplings one ideally needs samples of anomalous events, 
made by an event generator and one requires a fitting program containing the 
anomalous matrix element. The fitting program can then establish whether 
the input anomalous couplings can really be extracted from the generated 
anomalous data. 

Up to now such studies were made with tools, which have certain limi- 
tations. Usually data are generated for TT^-pair production containing three 
diagrams with PU-decay attached to it in zero width approximation. The 
fitting programs use the same approximation. Examples of such investiga- 
tions can be found in [13, 14]. Very recently a Monte Carlo program with 
anomalous couplings with a finite W- width became available [15]. It covers 
the semi-leptonic final states. 

In view of the advantages of the EXCALIBUR program it is natural to use 
its structure and strategy as a basis for an anomalous four-fermion generator. 
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Thus it is the aim of the present paper to describe the necessary additions 
and changes to the approach of [9, 10] to obtain a four-fermion generator 
with anomalous couplings. 

This generator has the following characteristics. Any massless four-fermion 
final state can be produced with CP conserving anomalous couplings. All 
abelian and non-abelian diagrams contributing to the final state are taken 
into account. It is also possible to restrict oneself to the "signal" diagrams of 
a process, e.g. for four-fermion production through IT^-pairs one takes only 
three diagrams. Finally, ISR can be switched on or off. 

The new anomalous matrix elements will be discussed in some detail, 
since they are the key ingredient of the anomalous EXCALIBUR program 
and since they would also be required for a fitting program. 

A number of numerical results for four-fermion production will be pre- 
sented. On one hand they serve as illustration how the W finite width or 
ISR can modify SM angular distributions. On the other hand they show how 
non-standard distributions behave. The very relevant physics application is 
the one mentioned above: generating anomalous data and studying the ex- 
traction of anomalous couplings with a fitting program. It is expected that 
in the future this question will be addressed. 

The actual outline of the paper is as follows. In section 2 the anoma- 
lous couplings are described. The next section discusses those four-fermion 
final states, which are sensitive to anomalous couplings and gives the re- 
quired matrix elements. Some illustrative examples of anomalous effects in 
distributions are shown in section 4, whereas section 5 contains conclusions. 

2 Anomalous couplings 

In this section those non-standard couplings are defined, which will be con- 
sidered for the generation of anomalous four-fermion final states. When one 
uses only Lorentz invariance as condition there exist 14 couplings, which lead 
to deviations from the SM triple gauge boson couplings. Some of them can 
be immediately discarded since they would either modify the strength of the 
electromagnetic interaction or introduce C or CP violation in it. At this 
point there still are 9 parameters left, three of which lead to CP violation 
through the ZWW interaction. Also these will be omitted. We are then left 
with a Lagrangian of the form: 

C = Ci(C- and P-conserving) + 

(1) 

£2 (CP-conserving, C- and P- violating) . 
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Here L\ takes the form 



(2) 



d = -ie [A^(W-^W + U - W^ V W- V ) + F^W+^W-"} 

+zecot 9 W [Z^W-^W + U - W + ^W-„) + Z^W^W-"} 

-iex^F^W+^W'" + iexzZ^W^W-" 

+ieS z [Z^{W-^W + V - W + ^W- V ) + Z^W + »W- V ] 

whereas the second part reads 



£ 2 = -||f d a Z pa (d p W- a W +a - d p W- a W +a + 

M w (3) 



dPW +a W~ a - d p W +a W~ a ) 



where Z^ is the dual field tensor 



Zpa — 7; e pcrapZ al3 (4) 



with Com = — 1- Taking the anomalous couplings vanishing leaves us with 
the first and second lines in Ci, i.e. the SM Lagrangian. Although this 
general form will be considered, there are theoretical symmetry arguments 
to reduce the number of independent couplings [13]. In practical fits this 
reduction will be necessary. For completeness we list the Feynman rule for 
the ZWW vertex when all particles are considered to be outgoing: 




= ie (cot w + 6 z )[g pv (p ~ P + )» + g»v(P + ~ P°) P + 9up(p° ~ P )•] 
+iex z [g^pPl - g^P° p ] 

[(PtPlPn - PVpPv) + (P° • P + )i9wPv - g P vP») 

+ (p° ■ p~){g P vPt - g^Pp) + (p~ ■ p + )(g^p° P - 9 w p° v )] 

1V1 W 

For the jWW vertex one has to replace cot 9 W + 5 Z by — 1, x z and y z by — x 7 
and — y 1 . The coupling z 7 is zero. 

With these vertices the matrix elements for four-fermion production will 
be evaluated. It should be noted that the form chosen in the interaction cor- 
responds to that of [13]. In the Z couplings the signs look different from [13]. 
This is however compensated by the vectorboson-fermion couplings, which 
differ between the two papers. For the SM we use here for the photon-electron 
vertex ie^ and for the Z-vertex ie^(v — cry 5 ) with a = —(4 sin^ cos 
and v = a(l — 4 sin 2 6 W ). In [13] the latter is the same but the photon vertex 
has opposite sign. 

3 The matrix elements 

In the literature [16] studies have been made of the effect of non-standard 
triple gauge boson couplings on the following gauge boson production pro- 
cesses 

(1) e+ e- -> W + W- , 

(2) e+ e- -> W ev e , 

(3) e + e~ — > Z v e v e . 

They are described by 3, 9 and 7 diagrams respectively of which 2, 2 and 1 
diagrams containing a triple gauge boson vertex. In practice these processes 
lead to four-fermion final states. For a specific four-fermion final state not 
only the "signal" diagrams of the above reactions contribute but also "back- 
ground" diagrams, of which some contain also triple gauge boson vertices. 
Thus the anomalous couplings can contribute to the background diagrams 
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as well. Tables 1-3 list the leptonic, semileptonic and hadronic final states 
which can originate from one of the above signals. Moreover, the number 
of abelian diagrams (N a ) and of non-abelian diagrams (N n ) is given. For 
the actual calculation of the matrix element we proceed as in [9]. We first 
repeat the SM calculation and then extend it to non-standard couplings. 
Although many diagrams can contribute to a specific final state there are 
only two topological structures (generic diagrams), given in fig. 1. In these 
diagrams all particles are considered to be outgoing. The actual Feynman 
diagrams will be obtained by crossing those electron and positron lines which 
were assigned to become the colliding e + e~ pair. In the Abelian diagrams 
the charges of the fermions determine the character of the two exchanged 
bosons, which may be W + , W~, Z or 7. In the non-abelian diagrams, two 
of the vector bosons are fixed to be W + and W~ , and the third one can be 
Z or 7. In this way we avoid double-counting of diagrams. In principle the 
particles and antiparticles can each be assigned in six ways to the external 
lines. For the non-abelian diagrams we get at most 8 diagrams and due to 
the specific final states considered we get at most 48 abelian diagrams. Fix- 
ing a specific four-fermion final state all possible assignments are tried. Only 
those which are allowed for by the couplings survive in first instance. Since 
also successively all helicity combinations are tried certain diagrams do not 
contribute as can be seen from the numerator of the generic abelian diagram: 



Here we have disregarded the particle/antiparticle distinction since it is al- 
ready implied by the assignment of the external momenta. The helicity labels 
A, p,a = ±1 determine the helicity of both external legs on a given fermion 
line. Using the Weyl-van der Waerden formalism for helicity amplitudes [17] 
(or, equivalently, the Dirac formalism of [18]), the expression A can easily be 
calculated. For instance, for A = p = a = 1 one finds 



A(+, +,+;!, 2, 3, 4, 5, 6) = 4(31)*(46) [(51)*(21) + (53)*(23)] (6) 



where the spinorial product is given, in terms of the momenta components, 



A(\,p,a;p u p 2 ,P3,P4,P5,P6) = 
= u x {pi)Yu\{p2) 

x u a (p 5 )Yu a (p 6 ) . 



(5) 



by 
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Figure 1: generic diagrams for four-fermion production. The fermion mo- 
menta and helicities, and the bosons are indicated. The bosons can be 
either Z, W ± , or 7; V can be either Z or 7. 

We denote the expression of equation 6 by A (l, 2, 3, 4, 5, 6). All helicity 
combinations can be expressed in terms of A , as follows: 

A(+ + +) = A (l, 2, 3, 4, 5, 6) A{ ) = A (l, 2, 3, 4, 5, 6)* 

A(- + +) = A (2, 1, 3, 4, 5, 6) A(+ - -) = A (2, 1, 3, 4, 5, 6)* 

A(+ + -) = A)(l, 2, 3, 4, 6, 5) A(- -+) = A (l, 2, 3, 4, 6, 5)* 

A(- + -) = A (2, 1, 3, 4, 6, 5) A(+ - +) = A)(2, 1, 3, 4, 6, 5)*. 

The numerator in the non-abelian diagrams can also be written in terms of 
the function A: 

U\{pi)-f a U X {p2) U p {p?i)lv.Up{pi) U a {p*>)7 v U a {ps) 

x { g ^( Pl + P2 y + <r (p 5 +P6)" + <r(p 3 +p 4 n 

= A(A,p,a;l,2,3,4,5,6)-,4(a,p,A;5,6,3,4,l,2) . (9) 

Thus, for massless fermions, every helicity amplitude consists of a sum of 
very systematic, and relatively compact, expressions. 

Extending now the jWW coupling with non-standard terms from the 
previous section we get as new numerator 
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[A(X, p, a; 1, 2, 3, 4, 5, 6) - A(a, p, A; 5, 6, 3, 4, 1, 2)] 

+x 7 B(A,p,ct;1,2,3,4,5,6) 

[(Pa + p 4 ) • (Pi + P2)B(\, a,p;l, 2, 5, 6, 3, 4) 
-(Ps + Pe) • (Pi + P 2 )5(A, p, a; 1, 2, 3, 4, 5, 6) 
+ (P3 + Pa) ■ (p 5 + Pe)£(p, A, a; 3, 4, 1, 2, 5, 6)] 

+^-C(A,p,a;l,2,3,4,5,6) 

and for the new ZWW -vertex 



-(cot 9 W + 5 Z ) [A(A, p, a; 1, 2, 3, 4, 5, 6) - p, A; 5, 6, 3, 4, 1, 2)] 
-x z B(X, p, a; 1, 2, 3, 4, 5, 6) 

-tSt [(Pa + P 4 ) • (Pi + P 2 )5(A, a, p; 1, 2, 5, 6, 3, 4) 
-(Ps + Pe) • (Pi + P2)5(A, p, a; 1, 2, 3, 4, 5, 6) 
+(P3+P 4 ) • (p 5 +Pe)B(p,X,a;3, 4,1,2,5, 6)] 

-^-C(A,p,a;l,2,3,4,5,6) 

+-gf£>(A,p,a;l,2,3,4,5,6) 

where the new functions B, C and D are defined as 

B(+, +, +; 1, 2, 3, 4, 5, 6) = 2(31)*(35)*(26)(34) 

+2(43)*(15)*(42)(46), 

C(+, +, +; 1, 2, 3, 4, 5, 6) = 2 [(13)*(23) + (14)*<24)] x 

[(35)*(45) + <36)*(46)] [(51)*(61) + (52)*(62)] , 
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£>(+,+,+; 1,2, 3, 4, 5, 6) = 

2z(35)*(35)*(34)(56) [(13)*(23) + (14)*{24)] 
-2z(46)(46)(34)*(56)* [(13)*(23) + (14)*(24)] (14) 
+2i(31)*{31)*{34){12) [(53)*(63) + (54)*(64)] 
-2z(42)(42)(34)*(12}* [(53)*(63) + (54)*(64)] . 



The expressions for B, C and D satisfy the same relations (8) as A. When 
equations (12)- (14) are denoted by B , C and D then 



B(+ + +) = 50(1, 2, 3, 4, 5, 6) B{ ) = B (l, % 3, 4, 5, 6)* 

B(- + +) = B (2, 1, 3, 4, 5, 6) B(+ - -) = B (2, 1, 3, 4, 5, 6)* 

B(+ + -) = B (l, 2, 3, 4, 6, 5) B(- - +) = B (l, 2, 3, 4, 6, 5) 4 

B(- + -) = B {2, 1, 3, 4, 6, 5) B(+ - +) = 5 (2, 1, 3, 4, 6, 5)* 



(15) 



and the same holds for C and D. 

Finally it should be noted that the vector boson propagators are imple- 
mented in the form (q 2 — My + iMyTy)~ 1 , irrespective whether q is timelike 
or not. This recipe guarantees the validity of electromagnetic gauge invari- 
ance. When this is violated even by a small amount forward electron cross 
sections can be off by orders of magnitude. This is due to photon exchange in 
the t-channel and was already noticed a long time ago as a problem in single 
VF-production [19]. A less ad hoc solution to this problem is underway [20]. 



4 Results 

Whereas at high energies total cross section measurements will give crucial in- 
formation on the size of possible non-standard couplings, one has to consider 
at LEP2 angular distributions for this purpose. The natural five-dimensional 
differential cross section is 

da(e-e+ — > W-\V + — > hhhh) (lg) 

d cos M cos ^id^id cos ^ 2 d^>2 

where 6 is the angle between the incoming electron and W~ . The angles 9i, 
cf>i are the polar and azimuthal angles of the particle f\ in the rest system 
of the parent particle W~ , whereas the angles 9 2 , 4>2 fulfill a similar role for 
the antiparticle originating from W + . The angles are defined with respect 
to coordinate frames related to the W~ and W + . The z-directions are the 
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directions of the momenta of the vectorbosons. The y-axes are defined by 
respectively p_ x q_ and q + xp + , where p_, p + denote the momenta of the 
incoming electron and positron and <f_, q + the momenta of the W~ and W + . 

In the zero width limit the above cross section is directly related to the 
helicity amplitudes for on-shell VF-pair production and functions describing 
the decay of the vectorbosons [12, 13]. In principle direct fits to the above 
cross section could be performed. In practice one- or two-dimensional dis- 
tributions will often be used. In the following we shall study da/dcos^, 
dcr/dcos#i, da/d^i. 

The main purpose of this section is to illustrate effects of certain phenom- 
ena which have sofar not been incorporated in anomalous coupling studies. 
These are the effects of the finite IT^-width, of ISR and of background dia- 
grams. It is useful to define a number of (differential) cross sections a eval- 
uated under different assumptions. In the first place we introduce SM cross 
sections a SM ,on, &sm, &sm,isr and a S M,aii which are respectively on-shell, 
off-shell signal cross sections (i.e. with three diagrams), the off-shell signal 
case with ISR and the cross section containing all diagrams. Furthermore, 
we define (Jan, &an,isr, an d (JAN,aii which are (differential) anomalous cross 
sections calculated with the three signal diagrams, without or with ISR, and 
with all diagrams without ISR. 

The following ratios give an illustration of the effects of the finite width, 
the ISR, background diagrams and of non-standard couplings: 

(17) 
(18) 
(19) 
(20) 
(21) 



Ri = 








R 2 = 


&SM,ISR 




&SM 


Rz = 


&SM,all 






R4 = 


can 


i 






i?5 = 


O 'AN, all 




&SM,all 



The reaction which we take as example is 

e~e + — > e~i> e iid. (22) 
The following input parameters are used 



a' 1 = 128.07 
sin 2 m = 0.23103 
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M z = 91.1888 GeV 
M w = 80.23 GeV 
T z = 2.4974 GeV 
rV = 2.08 GeV. 

where sin 2 9 W is chosen in such a way that combined with the above running 
a value the correct is obtained: 

G " = V2si7e w M^ (23) 

For the ISR the usual value of a is used. It should be noted that the 
above experimental values for the total widths are incorporated in the prop- 
agators. In EXCALIBUR the decay widths of the W into a lepton pair or 
quark pair are independent from the input total width. They follow from the 
other input parameters. Since one would like to have s-dependent widths in 
the s-channel and because this would violate gauge invariance the following 
practical procedure is used. The s-dependent widths can be transformed into 
a constant width [21]. When this constant width is used in both s- and t- 
channel gauge invariance is ensured in a simple way, which numerically agrees 
well with theoretically more sound methods [20]. Thus the calculations are 
performed with propagators (q 2 — M v + zMyfy) -1 , where 



My = M v /Jl + lv , (24) 



rv = rv/^/i + 7^ (25) 

V = Ty/My. (26) 

With these input values various differential cross sections have been eval- 
uated. The SM and anomalous cross sections with all diagrams have to 
be calculated with cuts avoiding thus the singularities due to the massless 
fermions. In order to make meaningful comparisons the cross section asu m 
R 3 has the same cuts. The imposed cuts are 



E e - tu4 > 20 GeV (27) 

|cos0 e - iM -| < 0.9 (28) 

|cosZ(mJ)| < 0.9 (29) 

m ud > 10 GeV (30) 

where 9 is the angle between the outgoing particles and the incoming beams. 
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In tables 4-7 total cross sections a an, &an,isr, ° an (cuts), o an, ail are 
listed for an energy of 190 GeV. The SM differential cross sections are given 
in figure 2. Differential cross section ratios are given in the form of histograms 
in figs 3-10. From R 1 it is seen that the inclusion of the finite width already 
changes the cos# distribution by a few percent. Comparing the on- and off- 
shell a an a similar angular modification arises [22]. Similarly the inclusion 
of ISR or background diagrams introduce even larger modifications of this 
angular distribution. In order to show the effects of the various anomalous 
couplings histograms of and R 5 are presented with values ±0.5 for every 
coupling successively, the others being zero at the same time. When doing 
the analysis with the three signal diagrams both for SM and non-standard 
couplings (i? 4 ) the overall picture is roughly the same as for the case where 
both cross sections contain all diagrams (R$). The effects of the anomalous 
couplings show up most clearly in the cos# distribution as can been seen 
when comparing to the pictures of the cos 9\ and 4>\ distributions. 

5 Conclusions 

With the extended EXCALIBUR program it becomes possible to study ef- 
fects of anomalous couplings in all four-fermion final states which receive 
contributions from non-abelian diagrams. In this way finite width effects of 
the vectorbosons are incorporated and studies of ISR and background dia- 
grams can be made. Up to 2 TeV the program works efficiently. For studies 
at higher energies the present phase space treatment of the multiperipheral 
massive vectorboson diagrams should be adjusted, which in principle does 
not pose any problem. For LEP2 this is not yet required. 

From the presented results it is clear that in particular the distribution 
in the W production angle 6 is affected by the finite W width, ISR and 
background. Also here anomalous couplings show up most clearly. The 
results of this paper give a quantative assesment of the above effects, which 
have hitherto not been considered in the literature. 
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label 


final state 


N a 


N n 


total 


signals 


1 


e + e~v e v e 


48 


8 


56 


1 2 3 


2 




14 


4 


18 


1 2 


3 


e~v e v T T + 










4 


u e e + fi~D,j, 










5 


is e e + r~i> T 










6 




17 


2 


19 


1 


7 


T+T~V T V T 










8 




7 


2 


9 


1 


9 












10 




17 


2 


19 


3 


11 


V e V e T + T~ 










12 


V e V e V e l> e 


32 


4 


36 


3 


13 




11 


1 


12 


3 


14 


V e V e V T V T 











Table 1: leptonic four-fermion final states in e + e collisions. 
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label 


final state 


N a 


N 


total 


signals 


1 


e~i> e ud 


16 


4 


20 


1 2 


2 


e~v e cs 










3 


v e e + du 










4 


v e e + sc 










5 


[i^v^ud 


8 


2 


10 


1 


6 












7 


[i^v^du 










8 












9 


T~v T ud 










10 


T~V T CS 










11 


r + v T du 










12 


T + V T SC 










13 


v e v e uu 


17 


2 


19 


3 


14 


v e v e cc 










15 


v e v e dd 


17 


2 


19 


3 


16 


v e v e ss 
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Table 2: semileptonic four-fermion final states in e + e collisions 



label 


final state 


N 


N n 


total 


signals 


1 


uudd 


33 


2 


35 


1 


2 


ccdd 










3 


udsc 


9 


2 


11 


1 


4 


dues 











Table 3: hadronic four-fermion final states in e + e collisions. 
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y/s = 190 GeV: 


a SM = 0.6490 ± 0.0010 pb 






only VFVF-diagrams, no ISR 














yz 


^ Zj 


-0.5 


0.6729 
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Table 4: oan succesively calculated for all anomalous couplings vanishing 
but one. Each second row gives the error to a an- 





y/s = 190 GeV: 


a SM = 0.5790 ±0.0011 pb 






only VFVF-diagrams, ISR 








S z 
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z z 
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Table 5: (Jan,isr f° r similar case as in table 4 
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190 GeV: 


a SM = 0.45970 ± 0.00097 pb 






only VFVF-diagrams, 


no ISR, cuts 
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Table 6: a an with the imposed cuts. Cases as in table 4 





= 190 GeV: 


a SM = 0.4705 ± 0.0010 pb 






all diagrams, no 


ISR, cuts 






Sz 
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0.5 


0.4824 


0.4651 
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0.4737 


0.4808 


0.4816 
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0.0010 


0.0010 


0.0010 
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Table 7: oan^u with the imposed cuts. Cases as in table 4 
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Figure 2: Distributions for the Standard Model 
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Figure 3: The ratios Ri,R 2 ,R3 
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Figure 4: The ratios Rz^Rz 
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Figure 5: The ratio 
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Figure 6: The ratio R 5 
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Figure 7: The ratio R4 
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Figure 8: The ratio R 5 
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Figure 9: The ratio R4 
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Figure 10: The ratio R 5 
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